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Abstract. One of J. Stuart Dowker's most significant achievements has been to 
observe that the theory of diffraction by wedges developed a century ago by Sommerfeld 
and others provided the key to solving two problems of great interest in general- 
relativistic quantum field theory during the last quarter of the twentieth century: the 
vacuum energy associated with an infinitely thin, straight cosmic string, and (after an 
interchange of time with a space coordinate) the apparent vacuum energy of empty 
space as viewed by an accelerating observer. In a sense the string problem is more 
elementary than the wedge, since Sommerfeld's technique was to relate the wedge 
problem to that of a conical manifold by the method of images. Indeed, Minkowski 
space, as well as all cone and wedge problems, are related by images to an infinitely 
sheeted master manifold, which we call Dowker space. We review the research in this 
area and exhibit in detail the vacuum expectation values of the energy density and 
pressure of a scalar field in Dowker space and the cone and wedge spaces that result 
from it. We point out that the (vanishing) vacuum energy of Minkowski space results, 
from the point of view of Dowker space, from the quantization of angular modes, in 
precisely the way that the Casimir energy of a toroidal closed universe results from the 
quantization of Fourier modes; we hope that this understanding dispels any lingering 
doubts about the reality of cosmological vacuum energy. 
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1. Introduction 

1.1. Basics 

A wedge is a region ft C R 3 bounded by two intersecting planes on which, say, the 
Dirichlet boundary condition is imposed. All the standard linear partial differential 
equations (wave, heat, Schrodinger, resolvent, . . . ) for all the standard fields (scalar, 
electromagnetic, spinor, . . . ) can be studied there, and their respective Green functions 
constructed, by standard but not totally trivial methods. 

As the prototype of a wedge problem we consider the conditions defining the 
"cylinder kernel", which provides [U EJ El H] the most direct way of calculating the 
vacuum expectation value of the energy of a quantized massless scalar field: 

T(t, x, x) is defined in R+ x x fi, (1) 

d 2 T 

— = -V 2 T for x in ft, (2) 

at 2 

T(t,x,x')=0 for x in <9ft, (3) 

T(0,x,x) = <5(x-x) for x in ft, (4) 

T(t, x, x') is bounded as t — » +oo. (5) 
Explicitly, in the three-dimensional case, we have 

ft = {(r, 9, z): < r < oo, < 9 < 9 , -oo < z < oo}, (6) 

„ d 2 I d I d d 2 

V = 1 1 1 , (7) 

6(x - x') = 1 5(r - r')5(9 - 9')5{z - z'), (8) 
r 

and an additional implied boundary condition, 

T(t, x, x') is bounded as z —> ±oo. (9) 
Often the two-dimensional reduction is of interest: 

ft = {( r , 9): < r < oo, < 9 < 9 }, (10) 

^2 d 2 1 d 1 d 

V = 1 1 , (11) 

dr 2 r dr r 2 d9 2 

5(x-x') = 1 5{r -r')5{9 -9'). (12) 
r 

Remarks: The corresponding Green function for all of R 3 in the role of ft is 

T ° {t ^ ) = n 2 (t 2 + \l-x>\ 2 ) 2 - (13) 
It is proportional to the t-derivative of the fundamental solution of the Laplacian in R 4 , 
with if = 0. It follows that the most relevant dimension is 4, the same as the space- 
time dimension of the quantum field theory, and we henceforth revise our terminology 
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accordingly. It is important to understand that in the cylinder-kernel formalism, t is 
"imaginary" time, related to the physical time x° by it = x° — x '. The term "cylinder 
kernel" refers to the domain R + x Q in R 4 , not to the cylindrical character of Q itself 
in the present work. 

The angle 9q is "good" if 9q — tt/N for some integer N. In other words, the wedge 
considered is one of 2N equal sectors into which Euclidean space is divided. (For mixed 
boundary conditions (Neumann on one side and Dirichlet on the other), 8 is good only 
if N is even (the number of sectors is a multiple of 4).) In such a case the problem 
can be solved immediately by the method of images: Replicate the singular source rt8J) 
by perpendicular reflection of x' through each boundary plane, with a minus sign for 
the Dirichlet condition but a plus for Neumann. Continue the construction through all 
2N sectors; the two directions of continuation meet consistently. The solution of the 
PDE in R 3 with this multiple source is easily constructed as a linear combination of 
copies of T with its source moved to each of the images. That function satisfies the 
desired boundary condition on each plane, so its restriction to fl is the desired Green 
function T. 

If the angle is "bad", the replicated wedges cannot match up to constitute Euclidean 
space. The constructed covering space is a cone, which is locally flat but has a singularity 
at the axis. The circumference of a circle centered at the axis is 9\ times the radius, 
where 9\ is 9q times the number of sectors in the construction. (The number of sectors 
in such a case is not uniquely determined. One could either take the simplest case, two 
sectors, or choose the value that makes 9\ closest to 2n.) To solve a wedge with a bad 
angle by images, therefore, one must know the cylinder kernel for the related problem 
on the cone. The latter satisfies the same equations @, Q, @, @ (fl now being the 
cone manifold), but the boundary condition ^ is replaced by the periodicity condition 

dT dT 
T(t, r, 9 + l5 z, x') = T{t, r, 9, z, x'), ^ (t, r,9 + 9 u z, x') = — (t, r, 9, z, x'). (14) 

Once this Green function has been found, it can be used in the role of T to construct 
the T for the wedge. 

In both cone and wedge, there is a technical issue about uniqueness of the 
solution, which is resolved by choosing the solution of minimal growth as the axis is 
approached [5j [6] . 



1.2. History 

At the end of the 19th century Sommerfeld EJ El HD] used the method of images to 
solve the problem of diffraction of waves by edges and wedges. (The seminal paper [7] has 
been translated with extensive commentary in [TT].) His and most other early work was 
concerned with the wave equation and with the solution corresponding to a particular 
incident wave, but of course the foregoing remarks about boundary conditions and 
images also apply there. At that time it was natural to think of the cone as a Riemann 
surface in the sense of complex analysis, rather than a locally flat Riemannian manifold 
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in the sense of differential geometry. Therefore, initially only values of Q\ of the form 
2ttM were considered, and hence wedges of angles 6 = ttM/N. Sommerfeld [7J was 
primarily concerned with a sharp edge (conducting half-plane), for which 8q = 2tt and 
hence a 2-sheeted Riemann surface (9\ = 4tt) divided into two sectors (N = 1) suffice. 
He gave mere indications of how to treat the general case 8q = irM/N, but his method 
was made explicit by Carslaw [12J, and the case of a right-angle wedge (6*0 = -jr, three 
sheets, four sectors) was worked out in detail by Reiche [THj. (The case 2n — 9 = |, 
M = 5, 2N = 6, is treated in the notes to [TT].) Sommerfeld [7J pp 346-347, 357] 
also showed how to construct solutions on an infinite-sheeted Riemann surface, but did 
not really explain how to use them to solve the original problem for # 7^ ttM/N (and 
the commentators in [11] describe this passage as "surrealistic"). He gave an explicit, 
though still very brief, prescription in P, p 38], as cited in Carslaw [H], and Wiegrefe 
[16] gave a more thorough treatment. 

Sommerfeld constructed solutions on Riemann surfaces by conformal mappings from 
the complex plane, but then he represented them by certain contour integrals, which 
are periodic with period Q\ . Carslaw [131 EH] noted that these integrals continue to 
satisfy the desired differential equations and periodicity even when 9i/tt is not rational. 
Knowing this, it is possible to dispense with the Riemann surfaces and obtain the 
solution of the wedge problem — for arbitrary 6 = |#i — from this periodic function by 
adding a single image term. Carslaw's papers [121 [13], [H] are a l so notable for emphasizing 
the generality of the method, beyond the wave equation. Other work on wedge problems 
during this era is reviewed in the introduction to [17] . which also constructs Green 
functions for the wave equation for arbitrary 9q in terms of standard special functions, 
without relying on either Riemann surfaces or contour integral representations. See also 
[18] [191 120] for applications of concepts of scattering and diffraction theory. 

Until 1977 the mathematically simpler problem of a cone was thought of by 
physicists merely as a means to the end of solving the problem of a wedge; it would 
never have occurred to Sommerfeld to speak of a "cosmic string", for example. But 
Stuart Dowker was well versed in this classic literature, and he realized that it could 
be instantly applied to problems that had arisen in quantum field theory in curved 
space-time [211 1221 123, I2H 122]- Again, the simplest problem was originally seen as 
a tool for studying something more recondite: In [2T1 [22] the emphasis is on the 
analytic continuation of a cone manifold to "Rindler space", the space-time visible to 
a uniformly accelerated observer; the periodicity is in the (imaginary) time coordinate 
and represents the effective temperature of the vacuum state. (In lieu of a long list 
of references, we refer to a previous review article [2S]-) Ten years later, however, the 
cone as a spatial manifold came to be regarded as physically realistic, as the idealized 
zero-radius limit of a cosmic string. (For subtleties in this interpretation see j27J 128]-) 
The Green functions and the vacuum expectation values of the stress tensor near a 
straight cosmic string were calculated in [231 [21] and many other papers, including 

In addition to moving the focus from wedges to cones, Dowker [21] introduced 
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a conceptual shift: The infinite-sheeted Riemann surface .Moo, previously regarded 
as a necessary evil required only for irrational cone angles, is now regarded as the 
fundamental manifold from which the others are built. (See also |4l] and references 
therein.) A ^-periodic image sum (of periodically displaced copies of a Green function 
on .Moo) yields the corresponding Green function on Ai 9l , the cone of defining angle 9i . 
This construction (which we discuss in more detail in section [4]) is precisely analogous to 
the creation of Green functions on an interval (or in a rectangular box) as periodic image 
sums of Green functions on the real line (or a higher- dimensional Euclidean space). A 
Green function for a wedge is then obtained by one more step of image summation, 
this time involving a single reflected copy (or, more generally, a sum over some even 
number of images, including the original source point). In the special case 6\ = 2tt, 
of course, the Green function for Minkowski space is recovered. Although the Green 
function on .Moo is harder to calculate than that on M.^ and perhaps not particularly 
easier than that on a general Aig 1 , the picture in which the covering space .Moo is the 
basic object is nicely consistent with contemporary notions of path integrals (as formal 
representations of exact solutions) and classical-path sums (as the basic framework for 
high-frequency and semiclassical approximations). It seems appropriate to refer to .Moo 
(including any extra spatial coordinates such as z in ([6])) as the Dowker manifold and 
to the corresponding space-time model as Dowker space. 

Dowker worked primarily with the same contour-integral representations as 
Sommerfeld and Carslaw. Our present exposition is strongly influenced by a later 
paper by Smith [30], which derives the cylinder kernels in total dimensions 3 and 4 
without overt use of such complex analysis ("overt" because Smith, and we, freely use 
special-function identities that appear in standard handbooks without inquiring how 
they were proved). Formulas in different dimensions are related in a way systematically 
investigated by Guimaraes and Linet [ID], which we shall explain in due course. Another 
method leading to the same results for the cone is that of Helliwell and Konkowski [29J , 
analogous to earlier work for the wedge included in the famous paper [42]. 

Finally, we must step back to an earlier paper by Lukosz [2] , which has been largely 
overlooked by authors in relativity until recently. This paper deals with the cylinder 
kernel (and vacuum energy) for a wedge. It is based on an insight similar to Carslaw's 
[13] [14] . Lukosz uses the method of images to solve the problem for the very special 
cases #o = k/N, as in section [i~T] He works the solution into a closed form (i.e., one 
no longer involving a sum over sectors) and then observes that the result satisfies the 
conditions of the problem for all values of 9$ . His method of course gives also the 
solution for an arbitrary cone. It can be checked to agree with those of [2"2"| I3TJ1 [3D] , all 
somewhat differently derived. 

Papers continue to be written about vacuum energy in the presence of either 
wedges or cosmic strings, the scenario being complicated or generalized in various 
ways of greater or lesser physical importance. An undoubtedly incomplete list is 
[131 SH US ESI 113 SHI SHI ISSl ED E2l [Ml [SH ISS|- The most experimentally relevant 
generalization is to the electromagnetic field, whose study goes back to Deutsch and 
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Candelas [12] and has been continued in (56J E3 EH1 EH] and numerous other references 
that can be traced back from (581 EH] ■ 

2. Formalism in cylindrical coordinates 

In order to calculate the vacuum expectation values of all the components of the stress- 
energy tensor, T^ v , we use a different cylinder kernel, T, which is an antiderivative of 
T as defined in section T satisfies equations (Jl]H([9j except that Q is replaced by 

^(0,x,x') = 5(x-x'). (15) 

In the calculations that follow, either "0" or an overdot is used to refer to components 
and derivatives with respect to the physical time coordinate. 

The stress-energy tensor formula for a massless real scalar field in flat space is [60] 

T, u = d^d v (j> - ^ Vtlv d x <f>d x <j> + £[r^<9 A <9 A (</> 2 ) - W0 2 )] 

= l -d^d v( t> - ^dfivt + /%^<9 A «9 A (0 2 ) - W0 2 )] (16) 

in Cartesian coordinates. Here £ is the usual curvature coupling parameter, but we 
write i = P + \ because then (3 = is the algebraically simplest case; thus (3 = — \ 
for minimal coupling and /3 = — ^ for conformal coupling in dimension 3 + 1. In the 



second version of ( 16 ) we have dropped a term that vanishes by virtue of the equation 



of motion, d\d x <p = 0. 

To calculate the components of the tensor in cylindrical coordinates, one can 



use (16) with rj and d reinterpreted as the metric tensor and covariant derivative. 
Alternatively, one can proceed as Schwartz-Perlov and Olum did for spherical 
coordinates [61)]: Without loss of generality, calculate at a point where the r, 6, and z 
unit vectors point along the x, y, and z axes. Define the "_L" components of tensors to 
be along the 9 direction, but with respect to an orthonormal basis, so that they have 
the same physical units as the other components; thus on any scalar function 

<9±0 = X d e <f>, (17) 
r 



and in (16) the metric is r/oo = — 1, i] rr = rj±± = i] zz = 1, all other components 0. The 
only remaining complication is that in the terms involving second-order derivatives the 
derivatives of the local basis vectors must be taken into account. That can be done 
by means of Christoffel symbols as usual in relativity or by direct calculation as in 
elementary vector calculus, with the results 

d ± 2 <j ) =^d r <j ) +^d e 2 ( j ) (18) 



and (31). Note also that 

de 2 {<p 2 ) = 2{de<pf + 2<pd e 2 ^ (19) 

so that any one of <9g 2 (0 2 ), (<9#0) 2 , and (j)dg(f) can be eliminated in favor of the other two. 
Continuing to follow [50], one then notes that, because of the symmetry of the situation, 
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the expectation value of 2 does not depend on time or (in the absence of end plates) 
the z coordinate, so do{<p 2 ) = d 2 {4> 2 ) = 0. For Dowker space or a cone (but not for a 
wedge) the expectation value of <fi 2 also does not depend on 6. All these simplifications 
yield 

t 00 = - #] - 2m-m4) + jW], (20) 

T rr = \[{d r( t>) 2 - 0<9 r 2 0] + 2/3[Vr0], (21) 
Zi T' 

T_l± = ^Q(0) - ^<R0 + 2/3[d r (<R0)], (22) 

t zz = \[{d z ct>) 2 - 4>d 2 z 4>] + 2t3[d r {<t>d r 4>) + -4>d r 4>], (23) 

2 r 

where 

Q{4>) = {de^f - 0«% 2 = 2(d e <P) 2 - \d e 2 {^). (24) 

Now we can rewrite the components of T^ u in terms of the cylinder kernel, T, rather 
than <fi. Since (O|0(x)0(a; / )|O) = — \T(x,x'), we have, for example, 

(0<9 r 0> = -\d r T, (09 r 2 0) = -|,9 2 T, ((9 r 0) 2 ) = -|9 r ^ r T, (25) 

where x' is set equal to x after differentiation. (Strictly speaking, we should symmetrize 
2d 2 as d 2 + <9 r 2 , etc., but it can be seen that this does not matter.) The derivatives for 
the other coordinates are computed similarly. The components of T^ v finally become 

(Too) = -\d 2 T + f3[d r d r ,T + dlT + l -d r T] (26) 
2 r 

(T rr ) =-\id r d r ,T-d 2 T}-^(3d r T, (27) 
(T ±± ) = ^d r T + [d e 2 T - d e d e ,T] - f3[d r d r ,T + 9 r T] (general) (28) 
= —d r T + ^dfT - (3[d r d r ,T + <9 r T] (symmetric case), (29) 

(T zz ) = -\[d z d z ,T - ,9 2 T] - (3{d r d r ,T + d 2 T + -d r T\. (30) 

Henceforth the expectation- value brackets around components of the stress tensor will be 
omitted for simplicity of notation; and later the notation T^ v will be abused once again 
after a vacuum subtraction. Again in these equations it is understood that formally the 
two space-time points are set equal after the differentiations. In reality these diagonal 
values are divergent, and either the points must be kept separated until after a vacuum 
subtraction has been performed, or some other kind of regularization must be adopted. 

In the interior of a wedge it is possible for T r ± to have a nonzero vacuum expectation 
value, but we shall not consider that component of the tensor further in this paper except 
to record the needed formula 

dj_d\\(f) = d\\d±(p = - d r de<p — \ (31) 

Here, to prevent ambiguity we write dn for the radial covariant derivative. (When it 
stands alone, d\\ is equivalent to d r : d\ 2 (fi = d 2 (fi, for instance.) 
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3. Calculation of T 

We now consider the problem of calculating cylinder kernels in polar coordinates. 
Various methods can be used. Guided by physical motivations, Smith [30] used one 
method in three (total) dimensions and a different one in four. We have carried out both 
calculations in both ways, and here we demonstrate how to do the four-dimensional case 
each way, to stress that they are mathematically interchangeable. Finally, we summarize 
the three-dimensional results and describe how adjacent dimensions are related. 



3.1. Homogeneous boundary-value problem; r -6 -z expansion 

The first method is to solve the boundary-value problem set up in section [T] (or its 
modification by ( 15 )). Thus T satisfies for t > the homogeneous 4-dimensional Laplace 
equation, 

d 2 T d 2 T 1 8T 1 d 2 T d 2 T 



+ 



+ 



+ 



+ 



0. 



dt 2 dr 2 r dr r 2 d9 2 dz 2 
along with a periodicity condition, T{9 + 9\) = T(9), and the intial condition 

T(0, r, r') = -S(r - r')5(9 - 9')5(z - z'). 
r 

Expanding T in a Fourier series in 9 yields 

oo 

T(t,r,9,z) = £ e me{ ^T„(t,r,z) 



(32) 



(33) 



(34) 



with 



T n {t,r, z) 



0i Jo 



-ine(|f) 



r ; T(t,r,fl,z) d9. 



Using these expressions in ( 32 ) , we see that 

d 2 T n d 2 T„ ldT n n 2 /2tt 



+ 



dr 2 



+ 



r dr r 2 \ 9\ 



T + 



d 2 T n 
dz 2 



From (33) and (35) we obtain 

1 fOi 



T n (0,r,z) 



0i Jo 
1 



d9e 



- me (ff )l$f r _ r ')s(e - 9')S(z - z') 
r 



—olr — r )o(z — z ). 
r 



(35) 



(36) 



(37) 



We now make a further step of variable separation, T sep (t,r) = T(t)R(r)Z(z), and 
we define A = 2 j 2L and occasionally suppress "n" in the notation. From (36) we get 



We let 



T" R" Z" IR! X 2 _ 
Y + ll + ~Z + rll~^~ ' 



T" Z" R" IR! X 2 



(38) 
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and therefore conclude that 



R" + 1 R' + (u 2 - % 



Z = e 



ikz 



R = 0, 

T = e , u = u + k . 



(40) 



(41) 



The appropriate solution of (40) is the Bessel function J\\\(ur) (and not Y\x\(ur) because 
of the minimal irregularity at r = 0). Now we can define the combined Bessel and Fourier 
transform 



oo roc 



ui't „ikz 



T n (t,r,z) = / udu / dkT(uj,k)J\\\(ujr)e w e 



When t = 0, from (37) we get 



e-^'6(r - r')5(z - z>) _ 
T n (0,r,z) = - — = P(r } z) } 

9\T 



and from (42) we see that 

T n (0,r,z)= I udu I dkT(u,k)J lM (ur)e tkz . 



OO POO 



(42) 



(43) 



(44) 



Therefore we can solve for T(u, k): 

_ 1 roc /• x 

T(u,k) 



2ll J -oo JO 



dz I r dr J\\\(ur)P{r\ z)e 



-ikz 



-,-ixe' 



'\ „—ikz' 



2n9 1 



From this one finally gets 
T(jt,r,9,z) = 
or, in full glory, 



J\x\{ur')e 



oo roo 



(45) 



u i 



du dk T n{Uik)J\ X \{ur)e ine e- wt e lkz (46) 



! °° roo roo 

T(t,r,9,z,r',9',z') = —— / uduj / dkJ\x\(ur)J\\\(ur') 

2n9 ln= _ OQ Ja J-oo 

x e iKe-e') e -u't e ik(z~z') _ 



(47) 



The formula for T is the same as for T except for a factor — u' in the denominator 
to implement the integration with respect to t: 

1 l ^~ > POO 

T(t,r,9,z,r',9',z') = - — £ e iX ^ / w do; J w {ur) J\x\(ujr') 

x dk(u 2 + k 2 )^ 2 e- tV ^ WI e lk{z - z '\ (48) 



By jnH (3.961.2)] the k integral equals 2K (u(), where ( = yr 2 + (z — z') 2 . Hence we 
have 



1 °° roo 

T(t,r,9,z,r',9',z') = -— £ e^'W w du; J, A |(u;r) ^(wrO^oK)- 

JO 



(49) 
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The integral over u can be computed [61, (6.522.3)], giving us the Fourier-series 
representation of T, 

T(t,r,9,z,r',9',z') = ~\ t e^-^f^Y"', (50) 

TtVi n i~' 00 r x r 2 \r 2 + r 1 J 



where r\ = J(r — r') 2 + ( 2 and r 2 = J (r + r') 2 + ( 2 . Let 

= e"«. (51) 

r 2 + r 1 

Smith [30] points out the alternative formulas for u that we present in section [4] and 
sums the series to get the final closed form of T, 

T(t,r,9,z,r',9',z') = - — 1 £ e^ 1 *"^ 

lixUirr sinn w 

1 sinh(|%) 



(52) 

27r#irr / sinhwcosh(§^) - cosg(# - 0') ' v ; 

(Actually, Smith calculated the Wightman function (a Green function for the wave 
equation) rather than T. It differs only by replacing t 2 by — (x° — x ') 2 .) 

3.2. Nonhomogeneous problem; t-0-z expansion 

The second method was employed by Smith [30J in dimension 3. In our case 
(dimension 4) it starts by observing that T (extended evenly in t) satisfies the 
nonhomogeneous Green-function equation 

d 2 T d 2 T ldT 1 8 2 T d 2 T 2 P/Nc/ u e/n „, u/ , ro . 
in all of R 4 . The proof of (53) is a standard application of Green's identity combined 



with an image construction in the t direction [62J. 



To solve (53) we again use a Fourier decomposition in 9, 



oo 



T{t,r,9,z) = E e m ^T n {i,r,z), (54) 



n=— oo 



T n (t, r,z) = i T e - me{ ^ } T(t, r, 61, z). (55) 
#i Jo 

Again let A = 2 W L ) and without loss of generality take 9=0 and z' = 0. By substituting 



into (53), we get 



d 2 T n d 2 T n ldT n A^ d 2 T n 2 , ,, . . , 

z _| z _| _ t H - = 5(t)5(r — r Wz). (56) 

dt 2 dr 2 r dr r 2 dz 2 9\r 

The delta functions in z and t can be represented as 

2 f 00 

8(z — z') = — dk cos kz cos fez', (57) 



7T JO 
2 /-oo 



8(t-f)=- du' cos uj't cos uj't'i (5£ 

7T JO 
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f n (u',r, k) 
T n {t,r, z) 



dt'cos(u't') / dz' ' cos(kz')T n (t' \r,z'), 



2 r°° 2 r°° ~ 

' du' cos(o/t) — / dk cos(fcz)T n (o/, r, k). 



TC JO 



TT JO 



Consequently, (56) becomes 



i '2 i 1 2\m i ® 2 Tn ldTn 2 , 

(w +^) T n + ^-^ + -^ ^T n =— 5(r-r). 

or z r or r z v\r 
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(59) 
(60) 

(61) 



Let or = lo + k (not the same relation as in (41)!). The solution of (61) has the 



form 



f 

J. , > 



C U I\\\ (or) for r < r', 
D^K^x^ur) for r > r'. 



(62) 



From the continuity and jump conditions implementing the delta function in (61), we 
get that 



= -jK\x\{ur') 



D u = -j-I\\\(ur') 



Therefore, 



and we get 



T n (u,r) = -— 7| A |(wr<)iif| A |(wr>), 



(63) 



(64) 



T(t, r, 9,z)= - — J] 



e e i 



oo 

OO POO 



x J dk J du' cos(co't) cos(kz)I^x\(cor < )K^x\(cor > ). (65) 

Since 2a/ da/ = 2a; da; and a;' = when u = k, the formula can be rewritten 

2 



e m( ^> uduI w (ur<)K w (ur>) 
i cos((o; 2 — k 2 )H) cos(kz) 

X / ClrC 



T(t,r,9,z) = 

,r(7 l n^oo JU 

■><■> (to 2 — k 2 )i ^ 

The k integral yields f JoOO, where ( = Vz 2 + 1 2 [61, (3.876.7)]. Thus by (HH 
(6.578.11)] 



T{t,r,9,z) 



^ r >)Jo{^0 



2tt\ fOO 



tt0i 



— v e ^(if) x e " l7r/2g w-i (coshM 



2tt rr'(smhu] 



1 5 
2 



(67) 



where 



coshw = 



r 2 + r' 2 + z 2 + t 2 
2rr' 
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and (5H (8.754.4)] 



i (cosh u) 



7T 



-\\\u 



2 sinh u 



Therefore, 



T{t,r,9,z) 



1ixQ\rr' sinhu n 



E 



^i\6-\\\u 



which agrees with the previous result (52). 



Remark: Several other variations are possible (see [63] for the 3-dimensional 
analogs). One can solve (53) by a complete Fourier analysis in all variables; then 



whether the first or second method above develops depends on whether the first integral 
one evaluates is over the variable u' conjugate to t or the variable u> conjugate to r. Yet 
another approach [121 !2H] is to introduce an additional "proper time" coordinate and 
find a Fourier solution for the corresponding heat kernel or propagator; then integrating 
over the proper time segues back to the approach just described, but evaluating the u' 
integral first yields a different representation. 



3.3. Three-dimensional problem; positive mass 

Both methods described above can be used in three dimensions, yielding the equivalent 
results [30] 



T(t,r,r',8) 



1 



E 



<2| A |_i(coshw ) 



71 



rr Ju 



du 



sinh(|^w 



-v/coshw — coshu cosh(|^w) — cos 



2tt6 



Here 



cosh M 



,2 + t 2 



A 



2nn 



2rr' 0i 
and Q is a second-kind Legendre function. For details see [30], [M] ■ 



(68) 
(69) 

(70) 



Formulas (68) and (69) cannot be simplified into an elementary closed form like 



(52) (just as a circular Bessel function is not elementary, but a spherical Bessel function 



is). However, (69) clearly is related to (52); indeed, the former is a certain integral 
transform of the latter. To understand fully the relations between cylinder kernels in 
different dimensions, it is necessary to enlarge the discussion to allow the scalar field to 
have a mass [10] ■ That is, a term m 2 T is added to the right side of In the case of 



two spatial dimensions, this would lead via (11) to the equation 
( d 2 



d 2 



d 



1 d 1 

dt 2 dr 2 r dr r 2 86 2 



m 



T = 



(71) 



(or the corresponding equation for T, or a similar nonhomogeneous equation like (53)) 



But if we merely rename m as k, (53) is recognized as the result of Fourier-transforming 



in z the four- dimensional massless field equation with spatial operator ([7]). Because we 



Wedges and Cones 



13 



are working in Euclidean signature when we study cylinder kernels, a massive Klein- 
Gordon equation is the same thing as a Helmholtz equation in one higher dimension. 

It follows that by taking a Fourier transform of a massive cylinder kernel with 
respect to the mass, one obtains a cylinder kernel for the massless theory in one higher 
dimension. (The massive higher-dimensional case can be reached by using only a piece 
of the lower-dimensional mass in the transform [401 (2.5)].) Conversely, by Fourier- 



transforming the four-dimensional T (52) in z, and then setting the Fourier variable 



equal to zero, one obtains the three-dimensional massless T; and (69) can be interpreted 



in this light after a change of variable from z to u (using (74), for example). 



4. The vacuum in periodic spaces 

In summary, we have (in total dimension 4) for the cone with angular periodicity 9\ 

1 sinh(^) 



T 01 (t,r, r',9,z) 



2tt9\tt i sinh u cos h 



2nu 



COS 



2nb 
01 



where u is defined by any of the equivalent formulas 

r 2 - r\ 



r 2 + T\ 



u — — hi 

r-i - 

T\ = \l{r — r') 2 + z 2 + t 2 , r 2 

/ i 2 /2 2 2 

2rr cosh u = r + r + z + t ; 



r + r') 2 + z 2 + t 2 : 



2rr' sinhw = J [r 2 + r' 2 + z 2 + t 2 ] 2 — 4r 2 r' 2 ; 



4rr' sinh 2 



u 



I /\2 i 2 i j.2 

r — r + z +t . 



(72) 

(73) 

(74) 
(75) 

(76) 



(Without loss of generality we have set z', and 9' equal to 0. When derivatives with 
respect to primed coordinates are needed, one simply replaces z by z — z', etc., first.) 



When 9 1 = 2tt, (72) reduces to the formula for Minkowski space in cylindrical 
coordinates, 

1 1 



T 



2tt 



An 2 rr' cosh u — cos 9 
1 



2vr 2 (r 2 + r' 2 + t 2 + z 2 - 2rr' cos(0)) 



(77) 
(78) 



(The two forms are shown equivalent by (74) and the law of cosines.) When 9\ — > oo 



(72) becomes the formula for the infinite-sheeted Dowker space-time, 

1 



T 



u 



(79) 

2ir 2 rr' sinhw u 2 + 9 2 

T$ 1 can be recovered from T m as a periodic image sum: after suppressing the 
irrelevant coordinates, 



9' + n9 1 ). 



(80) 
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In particular, 

oo 

7^(0 -00= E T oo (0-0 / + 27rn). (81) 

n=— oo 

In terms of normal modes, T2tt is a sum over angular momentum quantum number, while 
Tqo is an integral over angular momentum. T ro can be found by separation of variables 
(integration over modes) in a slight variation of the calculation in the previous section. 



Formula (80) can be verified, or discovered, in Mathematica or by complex analysis; 



indeed, the entire theory of vacuum stress in cones and wedges could be developed 



from formula (79) in this way. The Dowker Green function T ro is physically the most 
elementary; it tells how to diffract a path off the isolated conical singularity. The 
additional effects of periodicity, and hence of reflecting wedge walls, can be obtained 
from Too by images, or summation over paths, appropriate to each geometry in question. 
Of course, these remarks extend to the energy density and the rest of the stress tensor, 
and also to other Green functions, not just the cylinder kernel. 



The sum (80) has a structure similar to the formula for T in a universe that is 



periodic in one direction with period L± : 

oo 

T Ll (x-x')= ]T T (x-x' + nL 1 ). (82) 



Again, T^ is a Fourier sum, Tq a Fourier integral. The algebra leading from (79) to 



(72) according to (80) is very similar to that leading from (13) to a closed form for (82). 

There is, of course, a major glitch in the analogy of Too an d T dl to T and T Ll , 
respectively: it is T 27r (not Too) that, like T , is the free-space Green function that 
gives the zero-point energy density that must be subtracted from the energy of any 
other configuration. A periodic space has nontrivial Casimir energy, inherent in Tl 1 
after the subtraction. In the early stages of the universe, this vacuum energy should 
affect the cosmological expansion. 

Actually, this last conclusion is not universally accepted. The critique runs 
somewhat as follows: A real physicist is interested in the scalar field only as a model of 
the electromagnetic field. The (experimentally verified) electromagnetic Casimir energy 
is the energy of interaction of fluctuations of the electrons in the bounding conductors 
(the long-range limit of the van der Waals force). The field energy is just a bookkeeping 
device. But in a periodic universe (or a higher-dimensional closed universe) there 
are no boundaries, no fluctuating electrons. There is no experimental evidence that 
cosmological vacuum energy exists. (The observed dark energy may have something to 
do with zero-point energy, despite the notorious factor 10 120 , but it has nothing to do 
(in the absence of small extra dimensions, at least [65]) with the renormalized Casimir 
energy, which has a rest frame and is unobservably small at the present epoch.) 

Nevertheless, the mathematics (and indirectly confirmed physics) of quantum field 
theory unambiguously says that the vacuum energy density of a periodic universe differs 
from that of infinite space by a locally finite amount in precise analogy to how the 
vacuum energy density of infinite space (looked at in polar coordinates) differs from 
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that of Dowker space. In each case the increment is the difference between a Fourier 
series and a Fourier integral over the same spectral integrand. If you believe in the 
stress tensor of quantum field theory, this picture is totally consistent and unsurprising 
(except for the enduring mystery of what happens to the zero-point energy). But if you 
don't, you are forced into an untenable position: 

• To calculate vacuum energy in the periodic universe, you must ignore the 
(mathematically appropriate) Fourier sum in favor of the integral (since you don't 
believe vacuum energy can exist in the absence of van der Waals sources). 

• In polar coordinates, you must use the sum to get the right answer for empty 
Euclidean space. The integral gives something else, the energy density of the 
Dowker manifold. 

We see no possibility of a theoretical justification for this ad hoc switch of point of view. 
(Of course, this is not the same as an experimental or observational verification of the 
theory.) 



5. Numerics 



We now present some vacuum energy calculations for various cone and wedge spaces, 
based on the cylinder kernel for a cone as given in (72) and alternatively in (80), and 
executed in Mathematica. Most of these results were previously presented in pjo] . 

In this section we will work on diagonal, with r' — r, z — 0, and 9 = 0, but with 
t > 0. This is the traditional exponential ultraviolet cutoff. For flat space, we have 

Too= 2^> T - = T ^ = T - = ^I' ( 83 ) 
and of course we expect to subtract this "zero-point stress" from the stress tensor in 
any other configuration to get a physically meaningful quantity; this has been done in 
all the plots below, and in (84). Because t^T^ (with z = 0) depends on r and t only 
in the ratio r/t, the plot for t = 1 actually serves for all positive values of t, as the axis 
labels indicate. 

In the plots, solid curves are for the ultraviolet cutoff t ^ and dashed curves are 
for t = 0. When the cutoff is thus removed, the energy density and pressures are finite 
at each point in the interior of Q but (usually) develop nonintegrable singularities at the 
boundaries {9 — )■ 0, 9 — )■ 9 , or r — >■ 0). Note that when a cutoff is present (t > 0), the 
key quantity u is strictly positive even when all the coordinates are on diagonal (r = r', 
z = z' = 0, 9 = 9' = 0); therefore, all the quantities being plotted remain finite at the 
wedge plates, although the plots may give the impression of vertical asymptotes there. 

When multiple traces are shown in the same panel, blue, red, and (when present) 
yellow curves correspond to the cases in the order listed in the caption. To facilitate 
viewing on a monochrome device, we indicate in each case whether the corresponding 
ordering is from top to bottom or bottom to top in the bulk of the figure. Also, multiple 
plots within one figure are sometimes plotted with different horizontal and vertical scales, 
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in order to better show the behavior of the functions. In many cases, the energy density 
has been replotted in a separate graph to show the behavior for small r more clearly. 

Plots of the energy density and pressure in the Dowker space with /3 = (£ = |), 
obtained from (79), are given in figure [I] Figure [2] displays the additional terms in 
the energy density and pressure when /3 ^ 0. Because these correction terms are 
proportional to (3 = £ — j, it suffices to plot the case (3 = 1. 




a) Energy density 




(b) Radial pressure 




(c) Tangential pressure 



(d) Axial pressure 



Figure 1. Plot of the Dowker energy density and pressure as a function of r with 
t 



l 

4 ' 



In figures [3] and [| the energy density and pressure are presented for cone spaces 
with various cone angles. Notice that the overall sign of the energy density function 
changes at Q\ = 2tt. Plots for 9\ = 100007T match very well with plots for Dowker space, 
as expected since Dowker space is the large-period limit of the cones. The behavior of 
the energy density near r = is illustrated in figure |5j 

The curvature-coupling correction terms to the energy density and pressure are 
presented in figures [6] and [7] for various cone angles. The energy density for the conformal 
case, where £ = 1/6, is plotted in figure [8] and compared with the case £ = 1/4. 

In figures [9] and [TO] we present the energy density and pressure for the cone space, 
and the correction terms, as functions of the cone size Q\ with fixed r = 1. 
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(a) Energy density correction 




- 1 



(c) Tangential pressure correction 



(b) Radial pressure correction 




■ t 



(d) Axial pressure correction 



Figure 2. Plot of the Dowker energy density and pressure curvature-coupling 
corrections as a function of r with /3 = 1. 



For a Dirichlet wedge with opening angle 9 , we have an infinite sum over positive 
and negative images in Dowker space, and the resulting expression for T is 

1 / sinh(7ru/#o) 

47r#orr / sinh(w) \cosh(iiu/9 Q ) - cos(7r(6> - 6')/0 ) 

sinh(7ra/6'o) \ 



+ 



cosh(7ra/# ) - cos(tt(# + 6')/9 ) 
1 1 



4ir 2 rr' cosh u — cos 9 



34) 



The energy density now becomes a function of 9 as well as r. Figures 11 and [12] show 
the energy density as a function of 9 for various values of £, 9q , and r. Figures [lljb) and 
l~2~|b) show that when the cutoff is removed, the conformal energy density is independent 
of the angle (a fortiori, not divergent at the wedge plates), but at finite cutoff the function 
deviates strongly from this limiting behavior near the boundary. Both panels of figure 
11 show that the cutoff function is very far from the limit function when r is rather 



small. Figures [13] and [14] show that the energy density remains finite at the boundaries. 
Figures [15] and [16] show the energy density as a function of r for various values of £, 9 , 
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10M 



(b) Radial pressure 




(c) Tangential pressure 



(d) Axial pressure 



Figure 3. Energy density and pressure for £ = 1/4 and cone angles 0\ — 7r/4, 7r/2, 
and 7r (respectively blue, red, yellow, bottom to top in (a) and (c), and top to bottom 
in (b) and (d)). (Some plots for 6i = tt are missing because of a numerical problem.) 



and 6. 

Numerically, there is no qualitative difference between "good" and "bad" wedge 
angles, despite the great difference between them with regard to solvability by the 
classic method of images. The exception is the case 8q = it (for which the "cone" is flat 
space), where the "wedge" energy is just that of a single reflecting plane, expressed in 
polar coordinates. 

6. Conclusions and outlook 

There is a dual relationship between infinite (Euclidean) space and finite flat space (a 
circle or torus): Green functions for the infinite case can be obtained from the finite case 
by taking a limit, in which a Fourier series becomes a Fourier transform, or the finite, 
periodic case can be obtained from the infinite one as an infinite sum over image sources. 
A similar relationship exists in a two-dimensional space that is locally flat except at a 
single conical singularity (hence in a higher-dimensional space with a singular axis), 
except that the Fourier analysis is now with respect to the angular coordinate. The 
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r.'t ty't 



(a) Energy density (b) Radial pressure 




r.'t r/t 



(c) Tangential pressure (d) Axial pressure 



Figure 4. Energy density and pressure for £ = 1/4 and cone angles 9\ = 2.5tt, 8tt, 
and 10000-7T (respectively blue, red, yellow, bottom to top in (a) and (b), and top to 
bottom in (c) and (d)). 




(a) Energy density (b) Energy density 



Figure 5. Energy density for £ = 1/4 and cone angles (a) 6\ = tt/4, tt/2, and w 
(respectively blue, red, yellow, top to bottom), and (b) 9\ — 2.5ir, 8n, and 100007T (top 
to bottom). Behavior near r = 0. 
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01234: 01234? 
r.'t r,^ 



(a) Energy density correction (b) Radial pressure correction 




(c) Tangential pressure correction (d) Axial pressure correction 



Figure 6. Energy density and pressure correction terms for /3 = 1 and cone angles 
B\ = 7r/4, 7t/2, and ir (bottom to top). 

400 p — . — . — . — i — . — . — . — i — . — . — . — i — . — . — . — i — . — . — . — n 



300 - 
200 -. 

— 
£1 

100 - 

: L 




0.2 0.4 0.6 o.s l.o 

f : 

Figure 7. Energy density correction term for (3 — 1 and cone angles 9\ — 7r/4, tt/2, 
and 7r (top to bottom). Behavior near r = 0. 



role of the infinite space is now played by an infinite-sheeted singular surface, whose 
importance has been particularly stressed by Dowker. Among the continuum of periodic 
factor spaces (cones) of the Dowker manifold with respect to periodic image sums, one 
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o i i 3 4 s : •.: 02 OA 6 OS 1.0 

r/t r.i't 

(a) Energy density (b) Energy density 

Figure 8. Energy density for £ = 1/6 and £ = 1/4 (top to bottom in (a), and bottom 
to top in (b)), with cone angle 9\ = 0.8-7T. Behavior near r = is shown in (b). 




(a) Energy density 




(b) Radial pressure 




10 12 



(c) Tangential pressure 



(d) Axial pressure 



Figure 9. Energy density and pressure for the cone space as a function of 9i, with 
P = and r = 1. 



is our familiar, nonsingular Euclidean space. From the cone of arbitrary angle a further 
image sum, with respect to a single reflection, yields the Green functions for a perfectly 
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(a) Energy density correction 



(b) Radial pressure correction 




(c) Tangential pressure correction 



(d) Axial pressure correction 



Figure 10. Energy density and pressure curvature-coupling corrections for the cone 
space as a function of 9\ , with r = 1. 




(a) Energy density 



(b) Energy density 



Figure 11. Energy density in a wedge as a function of 9 with 9 = ir/2 and (a) 
£ = 1/4, for r at 2, 4, and 8 (top to bottom); (b) £ = 1/6, for r at 4, 8, and 16 (bottom 
to top). 
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(a) Energy density (b) Energy density 



Figure 12. Wedge energy density as a function of 6 with r = 8 and 6q values 7r/3, 
27r/5, and 2tt/3 for (a) £ = 1/4 (top to bottom); (b) £ = 1/6 (bottom to top). 




(a) Energy density (b) Energy density 



Figure 13. Wedge energy density as a function of 9 with 9 = ir/2 and (a) £ = 1/4, 
for r at 2, 4, and 8 (bottom to top); (b) £ = 1/6, for r at 4, 8, and 16 (bottom to top). 
Behavior near 6 = 0. 

reflecting wedge in Euclidean space, with an arbitrary opening angle. 

As shown by Lukosz and Smith, the particular Green function that is most useful 
for calculating the vacuum expectation value of the stress tensor in scalar quantum field 
theory can be found in closed form by rather elementary methods in all these spaces - 
Dowker space, cones, and wedges. The vacuum energy density (and pressure) in such 
systems consists of three components: the effect of the central singularity present already 
in the Dowker space, the effect of periodicity associated with a cone, and the specific 
effect of the reflecting edges of the wedge. The first two of these are independent of the 
angular coordinate inside the wedge. From this point of view, Euclidean space is itself 
a cone, whose periodicity energy precisely cancels the energy of the (erstwhile) central 
singularity. This periodicity energy is closely analogous to that of a torus relative to 
infinite space, which has sometimes been questioned because of the absence of sources. 
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Figure 14. Wedge energy density as a function of 6 with r = 8 and Oq values 7r/3, 
2ir/5, and 27r/3 for (a) £ = 1/4 (bottom to top); (b) £ = 1/6 (bottom to top). Behavior 
near 9 = 0. 




(a) Energy density (b) Energy density 



Figure 15. Wedge energy density as a function of r with 9o = tt/2 and 8 values 7r/16, 
7r/8, and ir/4 for (a) £ = 1/4 (top to bottom); (b) £ = 1/6 (bottom to top). 



We have presented numerous plots of the (renormalized) energy density and 
pressures in these systems, both with and without an exponential ultraviolet cutoff. 
In some circumstances a cutoff can be regarded as a mathematical device that can be 
completely removed at the end. However, alternatively a finite cutoff can be viewed 
as a small step toward a realistic physical model of the boundary material. In that 
case the standards for physical acceptability of the renormalized cutoff stress tensor are 
higher. Recent work on plane boundaries ^oT\ EE] indicates that the exponential cutoff 
introduces unphysical behavior in the energy density (but not the pressure) very near 
the boundary, and suggests that a better model of both energy and pressure could be 
obtained by replacing that cutoff (which corresponds to a separation of the arguments of 
the Green function in the imaginary time direction) by a point separation in a "neutral" 
direction. In the present cylindrical situation that would mean a separation in the axial 
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(a) Energy density (b) Energy density 



Figure 16. Wedge energy density as a function of r with 9 — it/2 and 8 values 7r/16, 
7r/8, and ir/4 for (a) £ = 1/4 (bottom to top); (b) £ = 1/6 (bottom to top). Behavior 
near r = 0. 



direction, and that prescription is currently under investigation, starting in [66J. 
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